First evidence of pep solar neutrinos by direct detection in Borexino 
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We observed, for the first time, solar neutrinos in tiie 1.0-1.5 MeV energy range. We measured the 
rate of pep solar neutrino interactions in Borexino to be (3.1±0.6stat±0.3syst) counts/(day-100ton) 
CO and provided a constraint on the CNO solar neutrino interaction rate of <7.9counts/(day-100ton) 

(95% C.L.). The absence of the solar neutrino signal is disfavored at 99.97% C.L., while the absence 
I of the pep signal is disfavored at 98% C.L. This unprecedented sensitivity was achieved by adopting 

novel data analysis techniques for the rejection of cosmogenic ^^C, the dominant background in 
y—^ the 1-2 MeV region. Assuming the MSW-LMA solution to solar neutrino oscillations, these values 

IL" correspond to solar neutrino fluxes of (1.6±0.3)x 10* cm"^s"^ and <7.7x 10** cm"^s"^ (95% C.L.), 

^ respectively, in agreement with the Standard Solar Model. These results represent the first mea- 

surement of the pep neutrino fiux and the strongest constraint of the CNO solar neutrino flux to 
date. 

PACS numbers: 13.35.Hb, 14.60. St, 26.65.+t, 95.55.Vj, 29.40. Mc 



Over the past 40 years solar neutrino experiments |T]- 
[5] have proven to be sensitive tools to test both astro- 
physical and elementary particle physics models. Solar 
neutrino detectors have demonstrated that stars are pow- 
ered by nuclear fusion reactions. Two distinct processes, 
the main pp fusion chain and the sub-dominant CNO 
cycle, are expected to produce solar neutrinos with dif- 
ferent energy spectra and fluxes. Until now only fluxes 
from the pp chain have been measured: ^Be, ®B, and, in- 
directly, pp. Experiments involving solar neutrinos and 
reactor anti- neutrinos [B] have shown that solar neutrinos 



undergo flavor oscillations. 

Results from solar neutrino experiments are consistent 
with the Mikheyev-Smirnov-Wolfenstein Large Mixing 
Angle (MSW-LMA) model [7], which predicts a transi- 
tion from vacuum-dominated to matter-enhanced oscilla- 
tioirs, resulting in an energy dependent Vf, survival prob- 
ability. Pee- Non-standard neutrino interaction mod- 
els formulate Pee curves that deviate significantly from 
MSW-LMA, particularly in the 1-4 MeV transition re- 
gion, see e.g. [8]. The mono-energetic 1.44 MeV pep neu- 
trinos, which belong to the pp chain and whose Stan- 
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dard Solar Model (SSM) predicted flux has one of the 
smallest uncertainties (1.2%) due to the solar luminosity 
constraint [9] , are an ideal probe to test these competing 
hypotheses. 

The detection of neutrinos resulting from the CNO 
cycle has important implications in astrophysics, as it 
would be the first direct evidence of the nuclear process 
that is believed to fuel massive stars (>1.5M0). Fur- 
thermore, its measurement may help to resolve the solar 
metallicity problem ^ [TU] . The energy spectrum of neu- 
trinos from the CNO cycle is the sum of three continuous 
spectra with end point energies of 1.19 (^'^N), 1.73 (^^O) 
and 1.74 MeV (^^F), close to the pep neutrino energy. 
The total CNO flux is similar to that of the pep neu- 
trinos but its predicted value is strongly dependent on 
the inputs to the solar modeling, being 40% higher in 
the High Metallicity (GS98) than in the Low Metalhcity 
(AGSS09) solar model [9]. 

Neutrinos interact through elastic scattering with elec- 
trons in the ^--^278 ton organic liquid scintillator target of 
Borexino [11]. The electron recoil energy spectrum from 
pep neutrino interactions in Borexino is a Compton-like 
shoulder with end point of 1.22 MeV. High light yield 
and unprecedentedly low background levels [51 [T^] make 
Borexino the only detector presently capable of perform- 
ing solar neutrino spectroscopy below 2 MeV. Its po- 
tential has already been demonstrated in the precision 
measurement of the 0.862 MeV ^Be solar neutrino flux 
[5J [131 ■ The detection of pep and CNO neutrinos is even 
more challenging, as their expected interaction rates are 
'^lO times lower, a few counts per day in a 100 ton target. 

We adopted novel analysis procedures to suppress the 
dominant background in the 1-2 MeV energy range, the 
cosmogenic /3+-emitter ^^C (lifetime: 29.4 min). ^^C 
is produced in the scintillator by cosmic muon interac- 
tions with ^■^C nuclei. The muon flux through Borex- 
ino is ~4300/^/day, yielding a ^^C production rate of 
~27counts/(day-100ton). In 95% of the cases at least 
one free neutron is spalled in the ^^C production pro- 
cess [13], and then captured in the scintillator with a 
mean time of 255 /iS [T5]. The ^^C background can be 
reduced by performing a space and time veto after co- 
incidences between signals from the muons and the cos- 
mogenic neutrons [TU [T7], discarding exposure that is 
more likely to contain ^^C due to the correlation be- 
tween the parent muon, the neutron and the subsequent 
"C decay (the Three-Fold Coincidence, TFC). The tech- 
nique relies on the reconstructed track of the muon and 
the reconstructed position of the neutron-capture 7-ray 
[T5] . The rejection criteria were chosen to obtain the 
optimal compromise between ^^C rejection and preser- 
vation of fiducial exposure, resulting in a ^^C rate of 
(2.5±0.3)counts/(dayl00ton), (9±1)% of the original 
rate, while preserving 48.5% of the initial exposure. The 
resulting spectrum (Fig. [l] top) corresponds to a fiducial 
exposure of 20409 ton-day, consisting of data collected 



— 10 
s 



Spectrum of events in FV 
Spectrum after TFC veto 

^^C rate = 27 ^^C rate = 2.5 

pep V rate = 3.1 CNO v limit = 7.9 

^^°Bi rate = 55 




0.8 1.0 1.2 1.4 1.6 

Energy / MeV 

FIG. 1. Top: energy spectra of the events in the FV be- 
fore and after the TFC veto is applied. The solid and dashed 
blue lines show the data and estimated ^'^C rate before any 
veto is applied. The solid black line shows the data after the 
procedure, in which the ^^C contribution (dashed) has been 
greatly suppressed. The next largest background, '^^"Bi, and 
the electron recoil spectra of the best estimate of the pep neu- 
trino rate and of the upper limit of the CNO neutrino rate are 
shown for reference. Rate values in the legend are quoted in 
units of counts / (day- 100 metric ton) . Bottom: residual energy 
spectrum after best-fit rates of all considered backgrounds are 
subtracted. The electron recoil spectrum from pep neutrinos 
at the best-fit rate is shown for comparison. 

between January 13, 2008 and May 9, 2010. 

The ^^C surviving the TFC veto is still a significant 
background. We exploited the pulse shape differences 
between e~ and e+ interactions in organic liquid scintilla- 
tors [TH1[I1] to discriminate ^^C /3+ decays from neutrino- 
induced e~ recoils and /3~ decays. A slight difference in 
the time distribution of the scintillation signal arises from 
the finite lifetime of ortho-positronium as well as from 
the presence of annihilation 7-rays, which present a dis- 
tributed, multi-site event topology and a larger average 
ionization density than electron interactions. An opti- 
mized pulse shape parameter was constructed using a 
boosted-decision-tree algorithm [5D] , trained with a TFC- 
selected set of -^^C events (e"*") and ^^'^Bi events (e~) se- 
lected by the fast ^^''Bi-^^^Po a-/3 decay sequence. 

We present results of an analysis based on a binned 
likelihood multivariate fit performed on the energy, pulse 
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Data (0.9 - 1.8 MeV) 




Pulse shape parameter 

FIG. 2. Experimental distribution of the pulse shape param- 
eter (black) . The best-fit distribution (black dashed) and the 
corresponding e~ (red) and (blue) contributions are also 
shown. 

shape, and spatial distributions of selected scintillation 
events whose reconstructed position is within the fiducial 
volume (FV), i.e. less than 2.8m from the detector center 
and with a vertical position relative to the detector center 
between -1.8 m and 2.2 m. We confirmed the accuracy of 
the modeling of the detector response function used in 
the fit by means of an extensive calibration campaign 
with a, P, 7 and neutron sources deployed within the 
active target |5J. 

The distribution of the pulse shape parameter (Fig. [2]) 
was a key element in the multivariate fit, where decays 
from cosmogenic ^^C (and ^°C) were considered e+ and 
all other species e^. 

The energy spectra and spatial distribution of the ex- 
ternal 7-ray backgrounds have been obtained from a full, 
Geant4-based Monte Carlo simulation, starting with the 
radioactive decays of contaminants in the detector pe- 
ripheral structure and propagating the particles into the 
active volume. We validated the simulation with calibra- 
tion data from a high-activity 228 ^j-^ source [21] deployed 
in the outermost buffer region, outside the active volume. 
The non-uniform radial distribution of the external back- 
ground was included in the multivariate fit and strongly 
constrained its contribution. Neutrino-induced e~ recoils 
and internal radioactive backgrounds were assumed to be 
uniformly distributed. Fig.[3]shows the radial component 
of the fit. 

We removed a events from the energy spectrum by 
the method of statistical subtraction [S]. We excluded 
from the fit all background species whose rates were es- 
timated to be less than 5% of the predicted rate from 
pep neutrinos in the energy region of interest. Further- 
more, we constrained all rates to positive values. The 
thirteen species left free in the fit were the internal ra- 
dioactive backgrounds ^lOfii, "C, ^"C, ''He, ^°K, ^SXr, 
and 234nipa (from ^ss^j jgcay chain), electron recoils 
from ^Be, pep, and CNO solar neutrinos, and external 
7-rays from ^ostI, ^^Bi, and ^or. We fixed the con- 
tribution from pp solar neutrinos to the SSM predicted 
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FIG. 3. Experimental distribution of the radial coordinate of 
the reconstructed position within the FV (black) . The best-fit 
distribution (black dashed) and the corresponding contribu- 
tions from bulk events (red) and external 7-rays (blue) are 
also shown. 

rate (assuming MSW-LMA with tan^ 6'i2=0.47l°;°4, 
Am^2=(7-6±0.2)xlO"'^eV2 i22J) and the contribution 
from neutrinos to the rate from the measured flux [J]. 
We fixed the rate of the radon daughter ^i^p]^ using the 
measured rate of ^i^'Bi-^^^Po delayed coincidence events. 

Simultaneously to the fit of events surviving the TFC 
veto, we also fit the energy spectrum of events rejected 
by the veto, corresponding to the remaining 51.5% of 
the exposure. We constrained the rate for every non- 
cosmogenic species to be the same in both data sets, since 
only cosmogenic isotopes are expected to be correlated 
with neutron production. 

Fits to simulated event distributions, including all 
species and variables considered for the data fit, returned 
results for the pep and CNO neutrino interaction rates 
that were unbiased and uncertainties that were consis- 
tent with frequentist statistics. These tests also yielded 
the distribution of best-fit likelihood values, from which 
we determined the p- value of our best-fit to the real data 
to be 0.3. Table |T] summarizes the results for the pep 
and CNO neutrino interaction rates. The absence of 
the solar neutrino signal was rejected at 99.97% C.L. us- 
ing a likelihood ratio test between the result when the 
pep and CNO neutrino interaction rates were fixed to 
zero and the best-fit result. Likewise, the absence of 
a pep neutrino signal was rejected at 98% C.L. Due to 
the similarity between the electron-recoil spectrum from 
CNO neutrinos and the spectral shape of ^lOgj^ whose 
rate is ^10 times greater, we can only provide an up- 
per limit on the CNO neutrino interaction rate. The 
95% C.L. limit reported in Table [l] has been obtained 
from a likelihood ratio test with the pep neutrino rate 
fixed to the SSM prediction [9] under the assumption 
of MSW-LMA, (2.80±0.04)counts/(day-100ton), which 
leads to the strongest test of the solar metallicity. For 
reference. Fig. |4] shows the full Ax^ profile for pep and 
CNO neutrino interaction rates. 

The estimated ^Be neutrino interaction rate is consis- 
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V Interaction rate Solar-!/ flux Data/SSM 

[counts/(day-100ton)] [10**cm"^s"^] ratio 
3.1 ±0.6etat± 0.3sy.t 1.6 ±0.3 1.1 ±0.2 

CNO < 7.9 (< 7.1,tatonly) < 7.7 < 1.5 



TABLE I. The best estimates for the pep and CNO solar neu- 
trino interaction rates. For the results in the last two columns 
both statistical and systematic uncertainties are considered. 
Total fluxes have been obtained assuming MSW-LMA and 
using the scattering cross-sections from |22H24j and a scintil- 
lator e" density of (3.307±0.003) x 10^^ ton"^ The last col- 
umn gives the ratio between our measurement and the High 
Metallicity (GS98) SSM [9]. 



Background 


Interaction rate 


Expected rate 




[counts/ (day- 100 ton)] [counts/ (day- 100 ton)] 


«^Kr 


19+^ 


30 ±6 [5] 


210Bi 








27.4 ±0.3 


28 ±5 


10^ 


0.6 ±0.2 


0.54 ±0.04 


«He 


< 2 


0.31 ±0.04 


40k 


< 0.4 






< 0.5 


0.57 ±0.05 


Ext. 7 


2.5 ±0.2 





TABLE II. The best estimates for the total rates of the back- 
ground species included in the flt. The statistical and system- 
atic uncertainties were added in quadrature. The expected 
rates for the cosmogenic isotopes ^^C, ^"C and ^He have been 
obtained following the methodology outlined in [25]. The 
expected ^'^'''"Pa rate was determined from the ^'^*Bi-'^^''Po 
measured coincidence rate, under the assumption of secular 
equilibrium. Ext. 7 includes the estimated contributions from 
20STI, ^^''Bi and *°K external 7-rays. 

tent with our measurement 5,. Table [TT| summarizes the 
estimates for the rates of the other background species. 
The higher rate of ^^°Bi compared to ,,5 is due to the ex- 
clusion of data from 2007, when the observed decay rate 
of ^^"Bi in the FV was smallest. 

Table |III| shows the relevant sources of systematic un- 
certainty. To evaluate the uncertainty associated with 
the fit methods we have performed fits changing the bin- 
ning of the energy spectra, the fit range and the energy 
bins for which the radial and pulse-shape parameter dis- 
tributions were fit. This has been done for energy spec- 
tra constructed from both the number of PMTs hit and 
the total collected charge in the event. Further system- 
atic checks have been carried out regarding the stability 
of the fit over different exposure periods, the spectral 
shape of the external 7-ray background and electron re- 
coils from CNO neutrinos, the fixing of ^^^Pb and and 
neutrinos to their expected values, and the exclusion 
of minor radioactive backgrounds (short-lived cosmogen- 
ics and decays from the ^"^^Th chain) from the fit. 

Table U also shows the solar neutrino fluxes inferred 
from our best estimates of the pep and CNO neutrino in- 
teraction rates, assuming the MSW-LMA solution, and 
the ratio of these values to the High Metallicity (GS98) 
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FIG. 4. Ax^ profile obtained from likelihood ratio tests be- 
tween flt results where the pep and CNO neutrino interaction 
rates are flxed to particular values (all other species are left 
free) and the best-flt result. 



Source 


[%] 


Fiducial exposure 


+0.6 
-1.1 


Energy response 


±4.1 


^^"Bi spectral shape 


+1.0 
-5.0 


Fit methods 


±5.7 


Inclusion of independent *^Kr estimate 


+3.9 
-0.0 


7-rays in pulse shape distributions 


±2.7 


Statistical uncertainties in pulse shape distributions 


±5 


Total systematic uncertainty 


±10 



TABLE III. Relevant sources of systematic uncertainty and 
their contribution in the measured pep neutrino interaction 
rate. These systematics increase the upper limit in the CNO 
neutrino interaction rate by 0.8 counts/(day-100 ton). 

SSM predictions [9]. Both results are consistent with the 
predicted High and Low Metallicity SSM fluxes assuming 
MSW-LMA. Under the assumption of no neutrino flavor 
oscillations, we would expect a pep neutrino interaction 
rate in Borexino of (4.47±0.05) counts/(day-100 ton); the 
observed interaction rate disfavors this hypothesis at 
97% C.L. If this discrepancy is due to oscillation to 
or we find Pee=0.62±0.17 at 1.44 MeV. This result is 
shown alongside other solar neutrino Pee measurements 
in Fig. [5] The MSW-LMA prediction is shown for com- 
parison. 

We have achieved the necessary sensitivity to provide, 
for the first time, evidence of the rare signal from pep 
neutrinos and to place the strongest constraint on the 
CNO neutrino flux to date. This has been made possible 
by the combination of the extremely low levels of intrinsic 
background in Borexino, and the implementation of novel 
background discrimination techniques. This result raises 
the prospect for higher precision measurements of pep 
and CNO neutrino interaction rates, if the next dominant 
background, ■^^"Bi, is further reduced by scintillator re- 
puriflcation. 

The Borexino program is made possible by funding 
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FIG. 5. Electron neutrino survival probability as a function 
of energy. The red line corresponds to the measurement pre- 
sented in this letter. The pp and ^Be measurements of Pee 
given in [5] are also shown. The measurements of Pee 
were obtained from [31 IH US] , as indicated in the legend. The 
MSW-LMA prediction band is the la range of the mixing 
parameters given in [22] . 

from INFN (Italy), NSF (USA), BMBF, DFG and MPG 
(Germany), NRC Kurchatov Institute (Russia), and 
MNiSW (Poland). We acknowledge the generous sup- 
port of the Gran Sasso National Laboratories (LNGS). 
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